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Q\ Abstract 
O 

Nuclear magnetic relaxation rate in iron-pnictide superconductors is calculated using the gap function obtained in a micro- 
{N| scopic calculation. Based on the obtained results, we discuss the issues such as the rapid decrease of 1/Ti just below the transition 
^ temperature and the difference between nodeless and nodal s-wave gap functions. We also investigate the effect of Coulomb inter- 
Q action on l/Ti in the random phase approximation and show its importance in interpreting the experimental results. 
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Iron-pnictide superconductors have been one of the central 
issues in condensed matter physics since its discovery. Many 
experimental and theoretical works have been carried out to elu- 
cidate the symmetry and the structure of the superconducting 
gap function since the information on the gap structure gives 
a vital clue to determine the origin of the superconductivity. 
Recently, it is recognized that the iron-pnictide superconduc- 
tors show varieties of the structure of the gap function! l'|, i- e -> 
nodeless behavior in some systems and nodal behavior in other 
systems. 

NMR measurement of 1/TiT is useful in determining such 
detailed structure of the gap function and many experimental 
results have been published^) [3] |4l [5] |6] . Thus, we calculate 
1/TiT based on some models for iron pnictides and compare 
the obtained results with experiments. Since many theoretical 
works exist[7 8, 9]|T0][TT|, we concentrate on the unresolved 
issues such as the rapid decrease of 1 jT\ just below the transi- 
tion temperature and the discussion on the difference between 
nodeless and nodal s-wave gap functions. 

For this purpose, we use the five-orbital Hubbard model pro- 
posed in ref IfTZl. Since we want to discuss the difference be- 
tween the nodeless s-wave and nodal s-wave, we use follow- 
ing two models. One is a simplified model downfolded from 
the LDA calculation on the LaFeAsO system. Here "simplify" 
means that we neglect the three dimensionality and take hop- 
ping integrals only up to the fifth nearest neighbors. Other is 
a three dimensional model of the LaFePO system. All results 
shown below are obtained with band filling n = 6.1. 

Assuming that the coupling between the conduction elec- 
trons and nuclear spins is diagonal in the orbital basis and using 
the constant form factor A(q) = 1, \jT{T can be written as 



Note that a and b denote the orbital indices. Without Coulomb 
interaction, we have 
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where q = (ia>„,q) and G a t,(q) {F a /,(q)) is the normal (anoma- 
lous) Green's function in the orbital representation. We also 
take acount of the Coulomb interaction within the random 
phase approximation (RPA), where the enhancement factor 
(1 - x(a)V s ) is included. Note that £{q) is the matrix form 
of eq. (iZb and V s is the spin vertex defined as V s . . = U (a = 
b — c — d), U' (a — c,b — d,a + b), Jh (a — b,c — d,a + c), 
/ (a — d,b = c,a + b), (others). In what follows, we use 
U = 1.0 eV , J H = 0.2 eV, U = U' + J H + J and J H = J. 

As is noted above, we describe the coupling between elec- 
trons and nuclear spins and the Green's function in the orbital 
basis. Thus, we describe the gap function also in the orbital 
basis. In this orbital representation, the gap function becomes 
matrix form as A(k). Limiting the discussion on the s-wave 
channel, physically important entries of A(£) can be written as 

&22(k)=ci2 + 2£>2(cos k x + cos ky) + 4c2 cos k x cos k v 

-45siak x siaky, (3) 
^33W=a2 + 2£>2(cos k x + cos ky) + 4c2 cos k x cos k y 

+46smk x sink y , (4) 

A44(£)=a4 + 2£>4(cos k x + cos k y ) + 4c4 cos k x cos k y , (5) 

A23(£)=2J23(cos k x - cos ky), (6) 

A24(£)=2/f>24(sin k x — sin k y ). (7) 

Note that subscript 2, 3 and 4 represent zx, yz and x 2 - y 2 or- 
bital respectively. To determine the parameters in eqs. ([3]|7|, 
we solve the linearized Eliashberg equation with the effective 
interaction calculated within RPA (we follow the formalisms in 
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Figure 1: T*T c /TiT without Coulomb interaction where T* is the T\ value at 
T c . Calculation is done with 192x192 (100x100x8) meshes in the momentum 
space for LaFeAsO (LaFePO) system. Upper (lower) panel shows the results 
without (with) Coulomb interaction. 



ref 02) )■ Then, we obtain the gap function with interband sign 
reversal witch can be parameterized as (02, bi, C2, S, 04, bn, 
c 4 , d 2 s, Via) = (0.0075, -0.038, -0.0155, 0.0075, 0.21, 0.11, - 
0.047, -0.0054, 0.045) for LaFeAsO and (0.060, -0.030, -0.01 1, 
0.0050, -0.072, 0.024, -0.0044, -0.022, -0.017) for LaFePO. 
The important difference between the two results is that A^k) 
has weaker k dependence in LaFePO and this weak dependence 
leads to the node when the gap function is written in the band 
representation! 12 1 . Namely, the gap function is nodeless in 
LaFeAsO and nodal in LaFePO. 

For studying the gap function below T c , we assume the usual 
temperature dependence, i.e., we use 



A(T, k) = ahk) tanh (l .74 y/T c /T - l) . 



(8) 



In this work, T c is fixed to a large value of 0.02 eV to avoid the 
numerical difficulty. Scaling factor a in eq. §8§ is determined 
to have 2A max /k B T c ~ 8 that is typical for iron-pnictide super- 
conductors with A max being the largest gap value on the Fermi 
surface. 

First, we show the results without Coulomb interaction in the 
upper panel of Fig. [TJ In this calculation, we use the damping 
factor y = 0.05 T c to discuss the clean limit (y is defined as 
ici) n — > a) + iy). Although the used gap functions have uncon- 
ventional property (interband sign reversal), the results show 
that the small coherence peaks appears in both of LaFeAsO and 
LaFePO system. (The coherence peak in LaFePO is neglegi- 
bly small.) Inset of the upper panel of Fig. [TJ shows that 1/T\ 
decreases slower in LaFePO system that has the nodal gap func- 
tion compared with LaFeAsO system as is expected. Actually, 
l/Ti in LaFePO decreases slower than T 3 just below T c . 

Next, we show the results with Coulomb interaction treated 



within RPA in the lower panel of Fig. [TJ The most important 
difference between results with and without Coulomb interac- 
tion is the disappearance of the coherence peaks. Here, we use 
the same damping factor y = 0.05T C as in the case of without 
Coulomb interaction, i.e., we assume the same strength of the 
impurity scattering. This means that we can eliminate the co- 
herence peaks without including the impurity effects. The rea- 
son for this disappearance is that Coulomb interaction makes 
the interband contribution to 1/T\ more dominant than the in- 
traband contribution and the dominance of the interband con- 
tribution leads to the disappearance of the coherence peak for 
the gap function with an interband sign reversal|[8] QTJ. An- 
other important point is that Coulomb interaction gives faster 
decrease of 1/T\ just below T c as can be seen by comparing 
the inset of upper and lower panel of Fig. [TJ Paying attention 
to the temperature just below T c , the result for LaFeAsO sys- 
tem shows more rapid decrease than T 6 , which is the exper- 
imintally claimed rapid decrease|4|, and the result for LaFePO 
system shows the decreasing rate close to T 3 . Note that such 
details depend on the value of a as is the case without Coulomb 
interaction. 

In summary, we have investigated the difference of the node- 
less and nodal s-wave gap function through calculating \/T\T 
using the microscopically determined gap functions. We have 
also demonstrated that Coulomb interaction can be a possible 
reason for the non existence of the coherence peaks in the clean 
limit. Further, we have shown the importance of Coulomb in- 
teraction in the discussion of the power of \/T\ just below T c . 
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